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Abstract 

Quantum mechanical models and practical calculations often rely on some ex- 
actly solvable models like the Coulomb and the harmonic oscillator potentials. 
The D dimensional generalized Coulomb potential contains these potentials 
as limiting cases, thus it establishes a continuous link between the Coulomb 
and harmonic oscillator potentials in various dimensions. We present results 
which are necessary for the utilization of this potential as a model and prac- 
tical reference problem for quantum mechanical calculations. We define a 
Hilbert space basis, the generalized Coulomb-Sturmian basis, and calculate 
the Green's operator on this basis and also present an SU(1,1) algebra asso- 
ciated with it. We formulate the problem for the one-dimensional case too, 
and point out that the complications arising due to the singularity of the one- 
dimensional Coulomb problem can be avoided with the use of the generalized 
Coulomb potential. 

PACS numbers: 03.65.Fd, 03.65.Ge, 03.65.Nk 
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I. INTRODUCTION 



Local potentials have been used to model the interactions of the subatomic world ever 
since the introduction of quantum mechanics. Some of these (like the Coulomb potential) 
do not differ essentially from the forces observed in nature, while most of them (like the 
harmonic oscillator, for example) represent approximations of the actual physical situation. 
The potential shape, defined by the potential type and the parameters in it is usually chosen 
in a way that reflects the physical picture our intuition associates with the problem; therefore 
we can define attractive or repulsive, short-range or long-range potentials, etc. 

Some of the potentials used in quantum mechanics are exactly solvable. This means 
that the energy eigenvalues, the bound-state wave functions and the scattering matrix can 
be determined in closed analytical form. The range of these potentials has been extended 
considerably in the recent years by investigations inspired by supersymmetric quantum me- 
chanics (SUSYQM), for example. (See Refs. |TJ for recent reviews.) Most of the solvable 
potentials are one- dimensional problems, which means that they are defined on the (—00, 00) 
domain (or on some finite interval of it), or they can be reduced to radial problems in higher 
spatial dimensions, in which case they are confined to the positive real axis. The simplest ex- 
actly solvable potentials are those belonging to the so-called shape-invariant class ||. Many 
of the most well-known potentials, like the Coulomb, harmonic oscillator, Morse, Poschl- 
Teller, etc. potentials belong to this class. Usually 12 potentials are mentioned to have the 
property of shape-invariance ||, although two of them are just alternative forms of other 
shape-invariant potentials. 

A rather more general potential type is that of the Natanzon potentials [f§] , the solutions 
of which can be reduced to a single hypergeometric or confluent hypergeometric function. 
These potentials are written as complicated six-parameter expressions of the coordinate, 
and their energy eigenvalues are given by an implicit formula. Although these problems are 
solvable exactly in principle, only some specific Natanzon potentials have been studied in 
detail 0-§J, because the formalism is getting more and more complicated with the increasing 
number of parameters handled. Some of the Natanzon potentials contain certain shape- 
invariant potentials as special cases: the Ginocchio potential || can be considered as a 
generalization of the Poschl-Teller potential, while the generalized Coulomb potential || is 
a simultaneous extension of the Coulomb and the harmonic oscillator potentials. There are 
also further potentials beyond the Natanzon class. These include, for example quasi-exactly 
solvable (QES) potentials ||, conditionally exactly solvable (CES) potential flO| and the 
SUSYQM partners of Natanzon-class potentials of which certain CES potentials form 
a subclass |T^| . 

The study of exactly solvable potentials represents an interesting field of mathematical 
physics in itself, but results from this area are also essential for the description of realistic 
physical problems. More and more interactions can be modeled by making advantage of 
the rather flexible potential shapes offered by exactly solvable potentials. The solutions 
of these can be applied directly, or they can be combined with numerical calculations. In 
the simplest case analytical calculations can aid numerical studies in areas where numerical 
techniques might not be safely controlled. This is the case, for example, when bound-state 
wave functions with arbitrary node numbers are required, for certain singular potentials, or 
for complex potentials []13f| . As the next level of complexity, analytical solutions can supply 
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a basis for numerical calculations. 

Virtually all quantum mechanical methods rely in some respect on analytically solvable 
potentials. Very often their wave function solutions are used as a Hilbert-space basis. More 
powerful methods can be constructed if we select a basis which allows the exact analytical 
calculation of the Green's operator of an analytically solvable potential. Some time ago by 
utilizing the Coulomb-Sturmian functions as basis functions the Coulomb Green's operator 
was calculated analytically and a quantum mechanical approximation method was developed 
||14| . Recently this method was extended for approaching the really challenging three-body 
Coulomb problem |T5[]. We note here that in this work the exact analytical representation 
of the Coulomb Green's operator on the Coulomb-Sturmian basis was essential. In our 
recent work |16| we showed that if on some basis representation the Hamiltonian possesses 
an infinite tridiagonal matrix (J-matrix) form then its resolvent can be given in terms of 
continued fractions which can be calculated analytically on the whole complex plane. This 
theorem was exemplified by the Coulomb and harmonic oscillator Green's operators. 

The aim of the present work is to demonstrate that various techniques developed for 
the simple shape-invariant potentials, like the Coulomb and harmonic oscillator potentials, 
can be generalized for certain Natanzon-class potentials too without making the formalism 
significantly more involved, and that these more general potentials can help our understand- 
ing of both practical and theoretical problems. The actual example we consider here is the 
generalized Coulomb potential |§, which is the member of the Natanzon confluent poten- 



tial class p7| . This potential is Coulomb-like asymptotically, while its short-range behavior 
depends on the parameters: it can be finite or singular as well at the origin. Its shape there- 
fore can approximate various realistic problems, such as nuclear potentials with relatively 
flat central part, or atomic potentials that incorporate the effect of inner closed shells by a 
phenomenological repulsive core. 

An advantageous feature of the generalized Coulomb problem is that the equivalent of 
the Coulomb-Sturmian basis can be defined for it, together with all its implications. For 
example, the method of calculating the matrix elements of the Green's operator as described 



in Ref. [16] for the Coulomb and the harmonic oscillator potentials can be extended to this 



case too, so that the above two problems appear in the present treatment as special cases. 

Using a particular differential realization, an SU(1,1) algebra can also be associated 
with this potential, which reduces to known SU(1,1) algebras in the Coulomb and harmonic 
oscillator limits. 

Another important feature of this potential is that it can be formulated as a one- 
dimensional problem, but it can avoid the pathologic features of the one-dimensional 
Coulomb problem, which arise due to the |x| _1 -like singularity. Since this potential can 
be chosen to be finite at x — 0, it offers a convenient basis for calculations which require the 
use of a one-dimensional Coulomb-like problem. 

The arrangement of this paper is the following. In Section [Tl] we give the basic formalism 
for the generalized Coulomb problem by recapitulating and extending former results. Sec- 
tion |TTT| contains the main results of the paper: i) the discussion of the Coulomb-oscillator 
connection using the two limiting cases of the potential, ii) the determination of the ma- 
trix elements of the Green's operator, in) the description of an SU(1,1) algebra associated 
with the problem and iv) discussion on the one- dimensional case, which elucidates the com- 
plications related to the singularity of the one- dimensional Coulomb problem. Finally, we 
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summarize the results and discuss further possibilities in Sect. 



II. THE GENERALIZED COULOMB POTENTIAL IN D DIMENSIONS 

Let us consider the Schrodinger equation in D spatial dimensions with a potential that 
depends only on the radial variable r. Separating the angular variables from the wave 
function we obtain the radial Schrodinger equation 

= (-^ + i(I + £^)(J + ^i) + «(r)) Vi(r) = c^(r) , (1) 

where u(r) = 2mh~ 2 V (r) and e = 2mh~ 2 E. The centrifugal term, which depends on the 
angular momentum Z originates from the kinetic term, i.e. from the .D-dimensional Laplace 
operator after separating the angular variables. The bound-state wave functions solving (|]) 
are normalized as 

/>rsn 

|^(r)| 2 dr = 1 . (2) 

Here we recapitulate the results obtained for the generalized Coulomb potential in Ref. 
by adapting it to the .D-dimensional case. Changing slightly the notation used in Ref. 
we define the generalized Coulomb potential as 

v(r) = + + + {p - - h c 



r 2V 2 /v 2 ' yr T yr T Ah(r)(h(r) + 6) 
q 3C 5C9 



fc(r)+0 16(/i(r) + 0) 2 16(/i(r) + #) 3 
where /i(r) is defined in terms of its inverse function 

r = r(/i) =C"3 ^tanh- 1 f^-A^'j + + fl))§ j . ( 4 ) 

The /i(r) function maps the [0, oo) half axis onto itself and can be approximated with 
h(r) ~ Car and /i(r) ~ Cr 2 /(48) in the r — > oo and r — > limits, respectively. 
Bound states are located || at 

e« = - jP„ , (5) 

where 

P. = ^((™ + /?/2) 2 + ^) -(n + ^/2)] . (6) 

The bound-state wave functions can be written in terms of associated Laguerre polynomials 

as 
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3+1 



Mr) = C^pT [ T{n + ^+l + p n 9) ) {h{r) + ^) i (Mr))^exp(-^Mr))Lr i )(p n / i (r)) 

(7) 

Potential (^) clearly carries angular momentum dependence: its first term merely com- 
pensates the centrifugal term arising from the kinetic term of the Hamiltonian. Its second 
term also has r _2 -like singularity (due to ft, -1 (r)), and as it will be demonstrated later, it 
cancels the angular momentum dependent term in the two important limiting cases that 
recover the D-dimensional Coulomb and the harmonic oscillator potentials. The third term 
of (Q) represents an asymptotically Coulomb-like interaction, while the remaining two terms 
behave like r~ 2 and r~ 3 for large values of r. 

In the remaining part of this Section we consider the specific case of D = 3 and I = 0. 
The S matrix of the generalized Coulomb potential can be derived in complete analogy with 
that of the Coulomb problem for D = 3. Although this can only be done exactly for / = 0, 
the singular term imitating the centrifugal term in Eq. (0) can be defined to be part of the 



potential. Using the notation of Ref. |T^] the equivalent of the general regular and irregular 
solutions are 

i f ( - G \k,r) = C-%(h(r) + 6) * (/i(r)) 2 ^ 1 exp(-~/i(r))$(~ +iu,P; ph(r)) (8) 

and 

f^ G \k,r) = p-fe^(/i(r) + 0)2(/ i (r))^ 1 exp(-|/i(r))^(^ + i^/?;p/i(r)) (9) 

where $(a, c; z) and ty(a,c;z) are the regular and the irregular confluent hypergeometric 
functions [|19j, and p and v have to be chosen as 

p = pW = -J|, * = = (io) 



The functions (H) and (||) reduce to (f^ c \k, r) and f^ c \k, r) of Ref. |IS| in the 9 — >• Coulomb 
limit. The S matrix for I = is then expressed as 

*<*) = (-d«+'^| . (id 

r(f -w) 

The extra phase factor appears because of the r _2 -type singular term, which is now defined 
to be part of the potential. In the Coulomb limit this expression becomes part of the 
centrifugal term, which is dealt with separately. 

The long-range behavior of potential (§) suggests its use in problems associated with the 
electrostatic field of some charge distribution. The deviation from the Coulomb potential 
close to the origin can be viewed as replacing the point-like charge with an extended charged 
object. The relevant charge density is readily obtained from the potential using 

" W = -8^ W ' (12) 
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Choosing D = 3 and neglecting the terms singular at the origin by making the choice I = 
and (3 = 1/2 or 3/2 we get 



h 2 2qC 20qC9-9C 2 81C 2 9 135C 2 9 2 

~ ~ 8nme^~ (h(r) + 9) 3 + 8(h(r) + 9) 4 + 16(h(r) + 9f ~ 32{h(r) + 9f 

2C^_ f h(r) V /2 / g 3C IhCB \ 

r \h(r) + 9J \(h(r) + 6) 2 8(h(r) + 9) 3 16(h(r) + 6)y l ' 1 ' 

In Figs. 1 and 2 we present examples for the actual shape of potential (|3[) and the corre- 
sponding charge distribution fll3|) for various values of the parameters. It can be seen that 
this potential is suitable for describing the Coulomb field of extended objects. It is a general 
feature of potential (|3]) that for small values of 9 a (finite) positive peak appears near the 
origin, which also manifests itself in a repulsive "hard core" , corresponding to a region with 
positive charge density (see Fig. 2). 



III. DISCUSSION 

A. The Coulomb and harmonic oscillator limits 

As described in Ref. ||, the introduction of the generalized Coulomb problem was in- 
spired by a simple method of transforming the Schrodinger equation into the differential 
equation of the confluent hypergeometric function. The bound-state solutions of these prob- 
lems contain associated Laguerre polynomials. The three most well-known examples of 
this type are the Coulomb, the harmonic oscillator and the Morse potentials, which are 
all members of the shape- invariant potential class ||. These three potentials are recovered 
by specific choices of the variable transformation reducing the Schrodinger equation to the 
differential equation of the associated Laguerre polynomials. This is also reflected by the 
structure of the function h(r) in Eq. (^): when h(r) is proportional to r and r 2 , one obtains 
the Coulomb problem and the harmonic oscillator potential, respectively. These limits can 
readily be realized by specific choices of the parameters in Eq. (Q): the first one follows 
from the 9 — > limit ||, while the second one is reached by taking 9 — ► oo, while keeping 



C/9 = C constant |20 



Besides taking the 9—^0 limit, the Coulomb problem in D-dimensions is recovered from 
Eq. (H) by the (3 = 21 + D — 1 and C"^q = 2mZe 2 /h 2 , choices: the third term of (||]) becomes 
the Coulomb term, the fifth one vanishes, while the other three all become proportional with 
r~ 2 and cancel out completely. 

In order to reach the oscillator limit one also has to redefine the potential (^) and the 
energy eigenvalues by adding q/9 to both. (This is equivalent to CD9 in the notation of 
Ref. [0.) This choice simply represents resetting the energy scale: e = corresponds to 
v(r — > oo) for the Coulomb problem, and to v(r = 0) for the harmonic oscillator. (Note 
that the energy eigenvalues also have different signs in the two cases.) Besides C/9 = C, 
the q = q/9 2 parameter also has to remain constant in the 9 — > oo transition here. The 
potential thus adapted to the harmonic oscillator limit reads 
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w ^ r) + < .._^( I + £zi)( I + £zl) + I) 

qhjr) 3C 1 5C 1 



4M0(i + ^) 

(14) 




The harmonic oscillator potential is recovered from (|i~4|) by the /3 = I + D/2 and Cg = 
(2muj/h) 2 choice. The two last terms in (14|) vanish, the first and the second cancel out, 
while the third one reproduces the harmonic oscillator potential. The new form of the energy 
eigenvalues is 

/ i \ 

/ / 1 n z \ 2 1 /5\ 

(15) 

which indeed, reduces to the e n = (2mu/H)(2n + 1 + D/2) oscillator spectrum in the 9 — > oo 
limit. The wave functions ([?]) are unchanged, except for the redefinition of the parameters. 

Note that the generalized Coulomb problem establishes a link between the Coulomb 
problem and the harmonic oscillator potential in different spatial dimensions. This is best 
seen by inspecting the wave functions (H). If the (5 — 1 parameter of the associated Laguerre 
polynomial is required to be the same in the two limits, we get an interrelation between the 
value of the angular momentum and the spatial dimension to be used for the Coulomb and 
the harmonic oscillator case: 

1° + ^- = 2l c + D C -I . (16) 



Considering 1° = 2l c , Eq. ([16| ) implies D° = 2D C — 2 which establishes link between the 
(D c , D°) =(2,2), (3,4), (4,6), (5,8), . . . etc. pairs. (Of these the (3,4) pair corresponds to the 



Kustaanheimo-Stiefel transformation |21|.) This case is called the "direct map" between the 
Coulombic and oscillator solutions in Ref. p2"f, while with the 1° = 2l c +X choice the "general 



map" can be recovered. These results suggest that the generalized Coulomb potential can 
be used to formulate a continuous transition between the Coulomb and harmonic oscillator 
potentials, as opposed to the usual procedure that employs a unique variable and parameter 
transformation to reach this goal. 



B. The generalized Coulomb Sturmian basis and the matrix elements of the Green's 

operator 

The generalized Coulomb-Sturmian equation which depends on n as a parameter, has 
similar structure to the eigenvalue equation in Eq. (ffl) with potential (0) 



3C 5C6 C{(3-\){(3 



XMn r) = \ — ^ — 1 — h ^ 

' [ dr 2 16{h{r) + 6) 2 16{h{r) + 6) 3 4h(r)(h(r) + 9) 
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and is solved by the generalized Coulomb-Sturmian (GCS) functions 



(r\n) =Mp,r) = (fj^J^) ' [ph{r) + pff)^{ph{r)f^ exp(-^/i(r))LM(p/i(r)) . 

(18) 

Here p is a parameter characterizing the generalized Coulomb-Sturmian basis. The GCS 
functions, being solutions of a Sturm-Liouville problem, have the property of being or- 
thonormal with respect to the weight function C^{h{r) + Introducing the notation 

(r\n) = (j) n (p,r)C^(h(r) + the orthogonality and completeness relation of the GCS 
functions can be expressed as 

(n\n) = 8 nn , (19) 

and 

oo oo 

1 = 2 \n)(n\ = \n)(n\ . (20) 

rc=0 n=0 

Straightforward calculation shows that both the overlap of two GCS functions and the 
(n'\Ho\n) matrix element can be expressed as a tridiagonal matrix, therefore the matrix 
elements of the e — H operator also have this feature: 

(n|e - H \n') = <W (^~ p ( 2n + P ~ p0) - ^ (-^ + (2n + (3) 

- <W+i (n{n + (3 - f + ^ ) - <W-i ((n + l)(n + /3))^ f + ^ 

\U2p 4 / \C2p 4 

This means that similarly to the D-dimensional Coulomb and harmonic oscillator potential 
the matrix elements of the Green's operator can be determined by using continued fractions, 
as described in Ref. ||16|| . The present results, therefore, extend the applicability of this 



method to a new potential problem. The formulae presented here reduce to those in [16 



in the appropriate limits discussed in Subsection limits. The role of the Coulomb-Sturm 



parameter b used in [16] is now played by C^p/2. 



C. Algebraic aspects 

Here we show that an SU(1,1) algebra can be associated with the generalized Coulomb 
problem. The generators of this algebra are 

- h + 9 - , f3. ? ? p. 

^-^^h-wm- (22) 
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an they satisfy the commutation relations 



[Ji, h] = -Us [4 Js\ = Ui [4 Ji] = ii 2 • (23) 

As can be seen from Eq. flT7|), J 3 is diagonal in the basis ( ]IBD with eigenvalues m = n + |. 
The elements of this basis can then be associated with the infinite-dimensional discrete 



unitary irreducible representation of SU(1,1) called ]23| discrete principal series Dj~, for 
which the allowed values of m are 



m 



-j,-j + l,-j + 2,... (24) 



with j being negative. It is natural then to identify j as j = — |. Direct calculations show 
that the ladder operators connect the neighboring members of this basis: 

J+MPi r) = (Ji + iJ 2 )0n(p, r) = [(n + l)(n + /3)] 1/2 n+1 (p, r) , (25) 

J_0 n (p, r) = (Ji - iJ 2 )0 n (p, r) = [n(n + /? - l)] 1/2 n _i(p, r) . (26) 
We find that the eigenvalues of the Casimir invariant 

are | (| — 1) = j(j + l), as expected, and that they set the strength of the fourth term in (|TTD. 
For 9 this is the only singular term and it behaves like ^r~ 2 with 7 = (/? — |)(/? — |) = 
4j(j + 1) + |. It is interesting to inspect the allowed values of 7 for the different unitary 
irreducible representations of SU(1,1). For the discrete principal series [3 > 1 holds, 
which always secures — | < 7, i.e. the potential has repulsive or "weakly attractive" [24 



-type singularity. For the supplementary series p3[ — | < j < holds, which results 



in — I < 7 < |. This is exactly the domain where both independent solutions are square 
integrable at the origin |24]| : for < 7 < | one of these vanishes at r = and the other 
one is infinite there, while for — | < 7 < for both solutions vanish at r = 0. ^From 
( [Of) it is seen that solutions regular and irregular at the origin correspond to (3 > \ and 
13 < This seems to indicate that one regular solution (with /3 > 1) is associated with 
while the second square integrable solution, which is either regular or infinite at the 
origin (depending on j3) might be related to the supplementary series, for which < (3 < 1 
holds. For the sake of completeness we note that for the continuous series [H C° k and C\J\ 



j = —\ + ik (k > 0, real) is valid, which results in the strongly singular 7 < — -7 case. 



The solutions then oscillate infinitely near the origin and are unbounded from below [24 
which can be interpreted as the falling of the particle into the center of attraction |25| . The 
situation is similar to that described in Ref. |26| for the v(x) = 7sin _2 x potential: the 
various unitary irreducible representations of the SU(1,1) spectrum generating group there 
also corresponded to different types of singularities. 

We note that we analyzed the singularities of the generalized Coulomb-Sturmian equa- 
tion (0) and its solutions (|T%D, but similar considerations of the physical potential (|3]) and 
its solutions can also be performed taking D = 3 and I = 0. The algebraic construction, 
however, does not apply to this latter problem. This is because the bound-state solutions 
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(0) pick up extra n- dependence through p n , which is not accessible for the ladder operators 
otherwise changing n as in Eqs. (^5]) and (|26|) . 

The present realization of the SU(1,1) algebra is a special case of that described in 



Ref. [17| in relation with the Natanzon confluent potentials. Considering the Coulomb and 



harmonic oscillator limits discussed in Subsection [111 A| and setting the dimension to D = 3, 



the generators reduce to the forms presented for the two problems separately in Ref. [27fl. 
We note that the spectrum generating algebra associated this way with the radial harmonic 
oscillator problem in three dimensions is different from the two-parameter realization of 



the SU(1,1) algebra discussed in Refs. PREUf, because the ladder operators there are linear 



differential operators and the Hamiltonian is related to the Casimir invariant, while here the 
ladder operators are second-order differential operators and the Hamiltonian is essentially a 
linear function of generator J3. 

Finally, we note that similar limiting cases of other Natanzon-class potentials can also 
be treated in terms of SU(1,1) algebras. In Ref. [^J an SU(1,1) algebra associated with the 
Ginocchio potential || was analyzed, and it was shown that in two different limiting cases, 
which both result in the Poschl-Teller potential, two different algebras can be recovered from 
the original one: the first being a spectrum generating algebra and the second a potential 
algebra. 



D. The one-dimensional case 

The formalism developed previously is valid for D = 1 too, nevertheless, some particular 
properties of one-dimensional problems justify a separate treatment of this case. First, the 
implicit definition of the h(r) function in Eq. (f|) has to be extended to negative values of 
r, which we now denote with x. Using the notation of (f4|), we can write that x = r(h) for 
x > and x = —r(h) for x < 0. The normalization of the wave functions in (0) also has 
to be modified with a factor of 2 -1 / 2 , accounting for the fact that the integration now runs 
from —00 to 00, but since the potential is symmetric, the bound-state wave functions are 
either even or odd, therefore |^>(:r)| 2 is an even function of x. 

For one- dimensional problems it is natural to set I to besides D — 1, which eliminates 
the centrifugal term in Furthermore, in order to avoid r _2 -like singularities at x = 
the second term in (|3|) also has to be canceled by setting (3 to either 1/2 or to 3/2. Ele- 
mentary calculations show that the latter choice corresponds to bound-state wave functions 
that vanish at x — 0, and essentially represent physical solutions of the problem in higher 
dimensions as well, while the former choice recovers solutions that do not vanish in general 
at x = 0. These two possibilities can naturally be interpreted as odd and even solutions of 
the one- dimensional potential problem. Furthermore, for x > the two types of wave func- 
tions can be rewritten into a common notation (up to a sign) by making use of the relation 
of associated Laguerre and Hermite polynomials, when the former ones have a = 1/2 or 



a = — 1 /2 as parameters |30 ] : 
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^ { N =1 \x) = ^ -(h(x) + 9^exp(-^h(x))H N ((p N h(x))^ 

2 n (r(M)r(f + i)(jv + 1 + ^) 5 2 

(28) 

For x < Othe bound-state wave functions satisfy 

* N ( } ~ \ -< =1) (*) for N = 2n + 1 (/? = |) . ^ 



In ([280 we defined pN as 




(30) 

which reduces to Pj-jvj, where even and odd values of iV have to be chosen with f3 = 1/2 and 

(3 = 3/2, respectively, and the integer part of N/2 corresponds to n used in p n in Eq. (|). 

An interesting aspect of this potential is that it remains finite at x = (v (0) = — q/9 + 
C/(86 2 )) for any finite value of 9, however, a narrow, finite peak appears in the 9 — > limit, 
which then becomes an attractive — 3/(16r 2 )-like singularity in the Coulomb limit. This is 
due to the fourth term in @ and it corresponds to a "weak" singularity in the sense that 
the center of attraction is not strong enough for the particle to become infinitely bound 
25| . This finite barrier arising for small, but finite 9 values also introduces the possibility 



of studying tunneling effects in symmetric potential wells. We also note that besides being 
finite at x = 0, potential (^) has continuous derivative there, as can directly be verified. 

Based on these features the D = 1 version of potential @ can be used to analyze 
the peculiarities of the one-dimensional Coulomb potential defined as = — e 2 /\x\. This 
problem has been the subject of intensive studies in the past couple of decades, but there is 
still some controversy in the interpretation of the results (see e.g. |H[ for a recent review). 
The unusual features attributed to this singular problem include degenerate eigenvalues K 



interpreted in terms of a hidden 0(2) symmetry ||33|| , an infinitely bound ground state |32 



and continuous bound-state spectrum J34]. The last two of these were later found to be 
based on unacceptable solutions of the Schrodinger equation |3^j3^|, while the unexpected 



degeneracy was explained by an impenetrable barrier at x = 0, which separates the problem 
into two disjoint, non-communicating systems with identical energy spectra |36| and makes 



even the concept of parity obsolete here Most authors discussing the one-dimensional 
Coulomb problem agree that the usual techniques of quantum mechanics alone in dealing 



with potentials are not sufficient in this case. In Refs. |38| for example self-adjoint extension 
of the relevant differential operator has been discussed. 

The one-dimensional version of potential @ can be chosen such that it becomes close 
to non-singular potentials used in the approximation of the true one-dimensional Coulomb 
potential. In fact, with appropriate choice of q and 9 any desired Coulomb asymptotics and 
v (x = 0) value can be generated. Fig. 1 shows potentials with rounded-off shape near x = 
{9 = 1, q = 0.5) and also ones close to the Coulomb potential with a cutoff — e 2 /(|x| + a) 
(9 = 1, q = 2.5). In contrast with these modified Coulomb potentials, all calculations 
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can be performed exactly with (^). In the 9 — > limit the generalized Coulomb potential 
recovers the one- dimensional Coulomb potential supplemented with the — jqX~ 2 term. This 
means that the one- dimensional Coulomb potential cannot be reached exactly, nevertheless, 
reasonable approximations of it can be given, as it will be demonstrated below. 

The odd solutions, of course, vanish at x — 0, while the even solutions have non-zero 
value there as long as 9 > holds. In the 9 — > limit ^>jv=2„ (0) varies with 1 / 4 , so the even 
solutions also tend to zero at x — 0. This is in accordance with the discussion in Subsection 
pi C| on weakly attractive 7r~ 2 type singular potentials on the half line: for — | < 7 < 
both independent solutions vanish at the origin, so the wave functions are necessarily zero 
at r = 0. If we try to extend the N = 2n solutions ( pS| ) in the 9 = Coulomb limit to the 
(—00, 0) domain we find that due to its x 1 / 4 type behavior at the origin the derivative of an 
even wave function would not be continuous anymore. 

It is also interesting to compare the Coulomb limit of the one-dimensional generalized 
Coulomb potential with the true one-dimensional Coulomb potential. Switching off gradually 
the r yr~ 2 term we find that the two types of solutions contain confluent hypergeometric 
functions of the type $(— n, c; z) with c = 1 + (1 + 47) 1 / 2 and c = 1 — (1 + 47) 1 / 2 . Starting 
with 7 = — we recover the c = § and | corresponding to the two c = (3 values assigned to 
the allowed solutions of the one-dimensional generalized Coulomb potential and its Coulomb 
limit. Taking the 7 — > limit results in c = 2 and c = 0. The first one of these corresponds 
to the allowed solutions of the one- dimensional Coulomb potential (see e.g. Ref. f39fl), while 
the latter one has to be rejected, because the confluent hypergeometric function $(a, c; z) is 
not defined for c = 0. This shows that even solutions which are allowed for the generalized 
Coulomb potential in one dimension do not have equivalents in the true one-dimensional 
Coulomb case. 

We also note that the true one- dimensional Coulomb potential can be reached in the 
9 — > limit setting /3 to 2 from the beginning. In this case, however, a repulsive 3/(16r 2 )- 
type singularity appears, which cancels out only for 9 = 0. This repulsive singularity then 
separates the problem into two parts (with x > and x < 0) from the beginning, so the 
potential is not a real one-dimensional one. The (3 = choice, in principle, also gives the 
same potential, however, the solutions associated with (3 = are not defined, as we have 
shown previously. Using (3 = e or 2 — e with e being small we can get as close to the true 
one-dimensional Coulomb potential as we wish, but the repulsive singularity is then present 
for all 9 values. 

The reflection and transmission coefficients can be analyzed using the asymptotic be- 
havior of the general solutions of the one-dimensional problem. These can be chosen to be 
even and odd functions of x. Making use of Eq. (^) the even and the odd solutions can be 
defined for x > in terms of cp( G \k, x) for x > 0, setting (3 = | and |, respectively; and their 
extension to x < can be given using a formula similar to Eq. (|29|) . The two solutions are 



interrelated by Eq. 6.3(3) of Ref. [0. Due to the symmetric nature of the one-dimensional 
potential (v(x) = v(—x)) it is enough to analyze the asymptotic behavior of the solutions 
for x —>■ 00: the x —>■ — 00 case follows naturally. Straightforward calculations show that the 
reflection coefficient is 

m) e-*»( rg + i») r(j + i„) \ 
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Strong reflection is found for potentials having a (finite) barrier in x = (like those in Fig. 
1 with 9 = 0.01), while more regular shapes (like that in Fig. 1 with 9 = 1 and q = 2.5, for 
example) give weak reflection. The k corresponding to maximal reflection increases as 9 get 
closer to (i.e. as the barrier at x = gets higher). |-R(fc)| is small for very large and very 
small values of k. Fig. 3 presents illustration for \R(k)\ 2 for various values of 9. Our findings 
on the reflection coefficient of the one- dimensional generalized Coulomb potential seem to 
support the existence of the space splitting effect [K7|] valid for the Coulomb potential on 
one dimension: in situations close to the true Coulomb case (e.g. 9 —>■ 0) \R(k)\ 2 is close to 
1, which corresponds to a nearly impenetrable barrier at x = 0. This is also confirmed by 
our analysis on other approximations of the one- dimensional Coulomb potential, i.e. those 
with (3 = e or 2 — e above. 

Besides its mathematical aspects the one-dimensional Coulomb potential has physical 
relevance too, in the description of the hydrogen atom in strong magnetic field |40|| , for 
example. In such practical calculations it is reasonable to use a non-singular Coulomb-like 
potential instead of the true one- dimensional Coulomb potential: the finite size of the nucleus 
can be a justification for this. This means that the one-dimensional Coulomb potential might 
not be sufficient in such calculations: the basis defined with it simply does not contain even- 
parity states. In practical calculations therefore the use of bases like that assigned to the 
generalized Coulomb potential is necessary. 



IV. SUMMARY AND OUTLOOK 

We analyzed the generalized Coulomb potential and extended results for the D- 
dimensional Coulomb and harmonic oscillator potentials to this potential too, which includes 
these two simple problems as special cases. The results presented here have both theoretical 
and practical significance. The theoretical aspects concern the treatment of Coulomb-like 
problems in one dimension. Due to the singular nature of the one- dimensional Coulomb po- 
tential the usual methods of defining a physically relevant basis are not sufficient in this case. 
The generalized Coulomb potential avoids these complications and it supplies a basis (in- 
cluding both even and odd functions) which can be used instead of the Coulomb-Sturmian 
basis, which is otherwise a rather helpful tool in higher dimensions, e.g. the three-body 



Coulomb problem has been approached this way |lq| . We also gave a unification of the 
SU(1,1) algebras associated with the Coulomb and harmonic oscillator problems separately. 

The practical aspects of our work include the introduction of the generalized Coulomb- 
Sturmian basis, which we used as a basis for analytic representation of the Green's operator. 
The generalized Coulomb potential (being the member of the rather general Natanzon po- 
tential class) offers potential shapes that can approximate various physical problems more 
accurately than the conventional textbook potentials. In particular, its long-range Coulomb- 
like behavior makes it a good candidate for describing the Coulomb field of extended (non- 
point-like) charge distributions. It can be applied either directly, or as a basis for numerical 
calculations. Its direct applications can be envisaged in describing muonic atoms, where 
it can model the Coulomb field of the nucleus. It can also be combined with short-range 
potentials in the potential separable expansion formalism (|TJ|,|T3J, for example. Then the 
Coulombic asymptotics of the interaction between the strongly interacting charged objects 
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could be taken care of by the Green's operator and only the short-range (nuclear) interaction 
would have to be approximated. 

Similar extensions of other relatively simple (shape-invariant) potentials also seem possi- 
ble. The example of the generalized Coulomb potential shows that the techniques developed 
for these well-known potentials can be extended to more general potentials without much 
technical difficulty The generalized potentials "inherit" many important features of the 
simple potentials, of which they are extensions, and, at the same time, they can be adapted 
to more complex physical situations, due to their more flexible shape. There are other 
examples confirming the expectations that seemingly complicated formulae can be evalu- 
ated analytically even for highly non-trivial potentials. In Ref. ||41|| , for example, potentials 
phase-equivalent with the Ginocchio potential || have been determined analytically, using 
the transformations of supersymmetric quantum mechanics. 

Based on the present results we believe that the Natanzon-class potentials, which have 
been analyzed almost exclusively in theoretical studies until now, can be used in practical 
applications too. 
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FIGURES 




Figure 2: The charge distributions corresponding to the potentials displayed in Fig. 1 for 0=1 and g=0.5, 1.25, 
2.5. (The largest q corresponds to the lowest curve.) 
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Figure 1: \R(k)\ 2 plotted as a function of 0. The notation x — k and y — 6 is used. For technical reasons the plot 
is made only from 6 = 10~ 5 instead of 0. 
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